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1. INTRODUCTION 

Variational methods have proven to be invaluable tools in theoretical physics and 
chemistry, both for bound state problems and for the study of collision phenomena. For 
collisional problems variational methods can be grouped into two types, those based on 
the Schrodinger equation and those based on the Lippmann-Schwinger equation. The 
Hulthen-Kohn 1 3 method belongs to the first type, and their modern development for 
electron-molecule scattering, incorporating complex boundary conditions, is reported in 
chapter 1 of this book by Rescigno et al. 4 An offshoot of the Hulthen-Kohn variational 
method is the variational R-matrix method. 5 6 In chapter 8 of this book Schneider 
presents a general discussion of the R-matrix method, including the variational R- 
matrix. 

The Schwinger variational (SV) method, which Schwinger introduced in his lectures 
at Harvard University and subsequently published in 1947,® belongs to the second 
category. The application of the SV method to e-molecule collisions and molecular 
photoionization has been reviewed previously. 9 12 The present chapter discusses the 
implementation of the SV method as applied to e-molecule collisions. Since this is not 
a review of cross section data, cross sections are presented only to serve as illustrative 
examples. 

In the SV method, the correct boundary condition is automatically incorporated 
through the use of the Green’s function. Thus SV calculations can employ basis func- 
tions with arbitrary boundary condition. This feature enables the use of an L 2 basis for 
scattering calculations, and provided the initial motivation for applying the SV method 
to atomic and molecular physics. 10 ' 13 The initial success led to the development of the 
iterative Schwinger method 14 which uses single-center expansion techniques and also an 
iterative procedure to improve the initial basis set. The iterative Schwinger method has 
been used extensively to study molecular photoionization. 15 For e-molecule collisions, 
it is used at the static-exchange level to study elastic scattering 16 and coupled with the 
distorted wave approximation to study electronically inelastic scattering. 

The Schwinger multichannel (SMC) method, originally formulated by Takatsuka 
and McKoy, 18 19 is the first modern computational method for e-molecule collisions ex- 
plicitly designed to treat the multicenter nature of a polyatomic target. When Gaussian 
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functions are used for the V basis and plane waves for the incoming electron, all two- 
electron repulsion integrals are calculated analytically. It is also a multichannel method 
and readily treats electronic inelastic scattering as well as polarization effects in elastic 
scattering. In addition, integrals involving the Green’s function, generally considered 
the bottleneck in the SV method, can now be computed efficiently and accurately us- 
ing the insertion method. 20 Alternatively, highly parallel computing can be effectively 
employed^ in the numerical integration over the final three- dimensions for this type of 
integrals. 12 The SMC method has been found to be robust, and has been employed to 
study molecules of various sizes and with different bonding types. For example, it has 
been used to study the electronic excitation of [1,1,1] propellane by electron impact, 21 
the largest polyatomic molecule studied so far by ah initio methods. It has also been 
used to study electron scattering from the ground state of BeCO, a molecule with a 
weak van der Waals bond, to simulate scattering from an adsorbate in a physisorbed 
system. 22 Since the SMC method is the most commonly used SV treatment for e- 
molecule scattering, the main topic of this chapter will deal with its implementation. 

Both the SMC and iterative Schwinger methods calculate fixed-nuclei, body-frame 
T-matrices. Domcke has demonstrated that nuclear dynamics in e-molecule collisions 
can be treated efficiently using projection operator formalism and Green’s function 
approach. Results of the fixed-nuclei SMC calculation can be coupled directly with 
Domcke’s treatment of nuclear dynamics. As an example, we present a calculation of 
vibrational excitations of N 2 by electron impact which combines these two treatments. 


2. THE LIPPMANN-SCHWINGER EQUATION AND THE SCHWINGER 
VARIATIONAL PRINCIPLE 


In the integral equation approach, we look for solutions of the Lippmann-Schwinger 
(LS) equation instead of the Schrodinger equation. Let the Hamiltonian of the electron 
-f molecule system be written as 



Here i,j and k y l sum over electrons and nuclei, respectively; Z * and m* denote the 
charge and mass of the kth nucleus. Also, atomic units will be used unless specified 
otherwise. Note that H is symmetric and its eigenfunction antisymmetric with respect 
to electron exchange. Rewriting H as 



H = H a + V, 

(2) 

with 

H 0 = Hm — ~ V^ +1 , 

(3) 

and 

* 1 n p 7 

v -y y' Zk 

fti l r w+i — Ti| |rjv + x-R t |* 

(4) 


The use of N + 1 to label the continuum electron is completely arbitrary. Also, H M is 
the molecular Hamiltonian with eigenfunctions <£ m , 


H^f^m — £ m $ m . ( 5 ) 
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The LS equation is given by, 


$(+) = S + G (+ V$ (+) . 


( 6 ) 


with 

lim {E -H 0 + i6)G w {E') = S(E - E'). (7) 

S — >0 

and S is the solution of the corresponding homogeneous equation. 1 / Here we use the 
boundary condition of an incoming plane wave and outgoing spherical waves, which 
has been incorporated into the Green’s function Due to this feature, it is not 

necessary to use trial functions which satisfy the correct boundary conditions in solving 
the LS equation. The Green’s function G ^ is expressed in terms of the molecular 
eigenfunctions and the interaction-free Green’s function. 


gM (e) = £ l*« 

m = l 


^ expfjl Smlrjy+i - tn+i '|) 
Ii'at+i — i*jv+i 1 


<*«!• 


( 8 ) 


where 

\kl + e m = E, (9) 

The summation m is over the entire molecular spectrum, including continuum states. 
When Em > E ) Ar m becomes imaginary, corresponding to a bound electron. Indeed, 
including the continuum spectrum of the target allows the outgoing electron to be 
different from the incoming electron, a necessary feature in order to account for electron 
exchange. 2 ^ 

Based on the LS equation for and the corresponding equation for with 

outgoing plane wave and incoming spherical waves boundary conditions, 

*(-) = s + gHvjH, 


the Transition matrix, or T-matrix, of a scattering process can be written as 

T mn = = (¥J\V - VG< + V|¥< + >). (10) 

It is readily seen that a combination of the above expressions 

T mn = (S m |V|*L + >) + (fWlVISn) - - FG( + )F|< + ))- (11) 

results in a T-matrix that is stationary with respect to first-order variations in 't (+) 

and respectively. 

6T mn = (S m |V|¥<, + > + *¥<+>> + (*i -) ff|Sn) - + ^L +) > - T -n 

= 0 , 

l / We choose to use a homogeneous solution which is the product of the target wave function and 
the free-particle wave function of the continuum electron, without antisymmetrization between the 
two. As pointed out in Ref. 24, it is not necessary to antisymmetrize both the initial and final wave 
functions explicitly. If ^(+> is antisymmetric, then the T-matrix element will automatically pick up 
the antisymmetric part of S. This convention is used here so that the same S will also apply to the 
SMC equation discussed in Sec. 3. 

2 i A different approach to account for exchange is to incorporate the Hamiltonian with permuted 
indices into the interaction potential, U = V 4- P{ jy + \(H — E). See Bransden et al. 
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and 


ST mn = + <*L _) + 69tf\v\ S n ) - \V - VG< + >V|¥<+>) - T mn 

= 0 . 

Equation (11) is referred to as the linear form of the Schwinger variational principle. 

In most e-molecule calculations, it is common practice to expand and 3^) in 
terms of a set of antisymmetrized N + 1-electron trial functions /, 

*i +) =EM*»)/r, (12) 

r 

and 

t’ = E^)/.. (13) 

» 

Eq. (11) gives 

T mn = J2 b r(Sr n \V\fr) + J2 C ’(f>\ V \Sr>)-Y l Y, b rC'(f.\V-VG^V\f r ). (14) 

r a r a 

The requirement that T mn be stationary with respect to first order variations in the 
expansion coefficients b r and c 3 , 


8T mn /db r = 0, 


dT mn /dc. = 0, 

gives the following relationships, 


b r = Y, D ''if’\V\Sn), 

a 

(15) 

C, = J2(Sm\V\fr)D r ., 

r 

(16) 

with 

(£"% = (f r \v -VG^v\f,). 

The resulting variational stable expression for T mn is 26 

(17) 

T mn = Y,{S m \V\f r )D r ,{f.\V\S n ). 

(18) 


ra 


Equation (18) has the advantage that it is independent of the normalization of 
and 

The advantages of the Schwinger variational method have been discussed prev- 
iously. 9 10 Because the correct boundary condition is already incorporated into the 
Green’s function and because the variational stable expression of the T-matrix in 
Eq. (18) is independent of normalization, it allows flexibility in the choice of a ba- 
sis. Furthermore, the wave function in Eqs. (11) and (18) always appears together with 
the potential, i.e., V|’l5 r £*‘)) or Thus a trial basis for needs only to cover 

the region of space where V does not vanish. This feature allows us to use an L 2 basis 
to represent 9^ when V is a short-range potential. 

Another advantage of the SV method is that the Schwinger variational principle is 
one rank higher than the Kohn variational principle. 27 28 If the same trial wave function 
is used, the Schwinger method should give a better converged result. It should be noted 
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that, in comparing the two methods, it is important to use the same trial basis because 
the Kohn method, which requires basis functions with the correct boundary condition, 
generally uses a trial basis different from the Schwinger method. 

A major drawback of the SV method lies in the difficulty in calculating the Green s 
function matrix elements. 9 The matrix element ($b' ) |yG (+) V|^^ +) ) in Eq. (11) or 
(/ r |VG (+) V|/,) in Eq. (17) are nine-dimensional integrals if V includes a two-electron 
potential. This obstacle has recently been removed when it is demonstrated that the 
insertion technique can be used with efficiency and good accuracy for the evaluation of 
these matrix elements. This will be discussed in Sec. 3.5. Alternatively, highly parallel 
computing can be employed for the 3-dimension numerical integration over the wave 
vector in these matrix elements. 

2.1. A Simple Example of Potential Scattering 

To illustrate the use of Eq. (18) and to demonstrate the convergence property of 
the SV method, consider the elastic scattering of an s-wave by a weak potential U. The 
initial state is described by the Riccati-Bessel function, 

S = jo(kr). 

Since the potential is weak, jo(kr) can also be used as a reasonable representation for 
the trial functions of $ + and Here the trial functions are of the form 

$(+) = b 0 jo(kr), 

vjd - ) = cojo(kr). 

Because the initial boundary condition is described by S in the LS equation, the coef- 
ficients b 0 and c 0 need not be set to 1. Instead, they are used as variational parameters 
and optimized. Note that the use of jo{kr) as a trial basis is an example of a basis 
which does not satisfy the correct boundary condition. Otherwise a combination of 
Riccati-Bessel and Neumann functions should be employed. Because $ (+) and , 4' ( 1 are 
expanded in a one term expansion, the inverse of the D-matrix, in Eq. (17), is just the 
inverse of the matrix element itself. From Eq. (18), the SV expression of the elastic 
X-matrix is given by 

_ (jo{kr)\U\jo{kr)){j 0 (kr)\U\jo(kr)) 

00 “ " (jo(kr)\V-UGMU\j 0 (kr)) 

Since U is a weak potential, we can write 

1 _ 1 (jo(Hlt/G(+>t/|jo(fcr)) , x 

(jo(kr)\U - {jo(kr)\U\j 0 (kr)) [ (j 0 {kr)\U \j 0 (kr)) 

Substituting the above into the expression for T 00 , we find 

Too = (jo(kr)\U\j 0 (kr)) + (. j 0 (kr)\UG^U\j 0 (kr )) + ... 

The first and second terms of Too correspond to the first and second Born terms, re- 
spectively. 

In Eq. (10) three non- variational expressions for T 00 are given. They give identical 
results if the exact and are used. Otherwise they give different values for the 
T-matrix. The first relation, 

Too = (5|C/|^ (+) ), 
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gives the Born approximation for 7'oo- 3 ' 


Too = (jo(kr)\U\j 0 (kr)). 

This result is inferior to the SV expression. Using the third relation in Eq. (10), 

Too = - U(7< + >U|*< + >), 

T 00 deviates even further from the variationally stable result, 

Too = (jo(kr)\U\j 0 (kr)) - (jo(kr)\UG M U\j 0 {kr)). 

Now the second Born term is subtracted from, instead of added to, the first Born term. 
The above exercise clearly illustrates the power of the variational method. 


3. THE SCHWINGER MULTICHANNEL METHOD (SMC) 

Currently, the most frequently used form of the Schwinger variational principle for 
e-molecule collisions originates from the work of Takatsuka and McKoy. 18 19 As seen be- 
low, while this formulation is based on the Schwinger variational principle, its equation 
of motion is obtained by combining the LS and Schrodinger equations. These equations 
are coupled via the introduction of a projection operator P which projects into the open- 
channel target space. In the following, we shall first consider some of the properties of 
Takatsuka s projection operator, then derive the SMC equation. Implementation of the 
method will then be delineated. Efficient evaluation of the Green , s function matrix el- 
ements, generally considered a bottleneck in the Schwinger method, will be discussed 
and extension to include correlated target functions will also be presented. 


3.1. The 7V-Electron Projection Operator P 

The projection operator P introduced by Takatsuka and McKoy projects into the 
open channel target space. Hence it is an TV-electron projection operator and different 
from the Feshbach projection operator, 29 which is a N -f 1-electron operator. 

M 

P= £ |<M1,2,...,A))(<M1>2,...,IV)| I (19) 

m=l 


with m summing over all energetically accessible target states. P satisfies the idempo- 
tent property, 

P 2 = P, 


and commutes with H 0 


PH 0 = H 0 P. 

To illustrate the difference between P and P/v+ i ? the Feshbach projection operator, 
consider the simplest case when the target wave function is expressible in terms of 
a closed shell, single determinental wave function, 

4>j(l,2, • - - , N) = A ^4>iGc(l)4>if3(2 ) . . . 4> N j 2 a{N — 1)^jv/2/3(-AT)J. 

3y/ Here 'jd+l is not optimized and 6 0 = 1 is used. 
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Here a and j3 denote spin functions and A the antisymmetrizer. If the wave function 
of the continuum electron, g t {N + 1), is orthogonal to all the target orbitals, then the 
S z — 1/2 component of becomes 

$i(l,2 ,...,N,N + 1) = A^4>xa(\)4>^{2) . . .<j>Ni2p{N) gi a(N + l)j. 

We have 

P*i(l, 2 ,...,N,N + l) = ^====*i(l,2 ,...,N)g,a{N + 1), (20) 

and ^ 

(9i\BP\9i) = (21) 

On the other hand, the Feshbach projection operator is defined in terms of the N + 1- 
electron open channel functions, 

Pjv+1 = £ |» m (l, 2, . . . , N, N + l))<« m (l, 2, ...,N,N + l)\. 

m 

where satisfies the correct asymptotic boundary condition. The function tf; is an 
open channel function in the terminology of the Feshbach projector formalism if g, 
satisfies the correct boundary condition. In that case, 

Pn+1% = ( 22 ) 

A comparison of Eqs. (20) and (22) shows how the two projection operators differ. Also 
note that Takatsuka’s projection operator treats S'; as an open channel function as 
long as it is associated with an energetically accessible target function, regardless of the 
behavior of g, at the boundary. The Feshbach projection operator, on the other hand, 
includes 'I', in the open channel space only if g t has the correct boundary behavior. 

It should be noted that Eq. (20) may not be applicable when the orthogonality 
constraint between the continuum electron function and the target function is relaxed, 
and/or when the target is described by more sophisticated wave functions. Consider 
the simple case of X J E+ -> fc 3 E+ excitation of H 2 by electron impact. Let the wave 
function for the X‘E+ state be represented by 

= ^{lc 9 a(l)l<Tj/?(2)}, 

and the three spin components of the state by 

< I , 6 1 s I =i = • / 4{l (7 9 a (l)l <T « a (2)}> 

^6,s,=o = «4{1 o- 9 (1)1<Tu( 2)— ;= [a(l)/3(2) + /3(l)a(2)] }, 

= A{la g /3( l)l ffu /?(2)}. 

It is well established that in calculating this transition, terms of the form 

= ,4{l<7 9 a(l)l<7 u a(2)l<7„/?(3)}, (23) 

and 

= A{ l<7 9 a(l)l<7 9 /3(2)l<7 u a(3)}, (24) 
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should be included in the expansion of the wave function. These are usually called ‘pen- 
etration’ or ‘recorrelation terms’ because they relax the enforced orthogonality between 
the continuum and target functions. Let the projection operator P be generated from 
the four target functions above. We then have 

= ^M.=i(1,2)1<t u /3(3) - -J=$ t ,s, = o(l,2)l<7 u a(3). (25) 

This simple example illustrates how the operation of P depends on the tvpe of functions 
used. 


3.2. The SMC Equation 

The projected LS equation is 

^l +) = S n + <#V*<+>, (26) 

The projectei. Green’s function, G ( p\ is defined in the open channel space by, 


gJ, +) 


2tt , /i 


\ r N+i — r N+1 'j 

Multiplying Eq. (26) from the left by V and rearranging, we find 

(VP - VG<p + W)*M = VS n (27) 

Equation (27) describes only the open channel functions, whereas a complete de- 
scription of ¥<+> also requires the closed channel component. To do this, we use the 
following identity 

a.P +(l - aP) = 1. ( 28 ) 

Note that in association with the operation of P, a parameter a, called the projection 
parameter, is introduced. It is seen from Eqs. (20) and (25) that the operation of P 
not only removes the antisymmetrization between the target function and continuum 
orbital, but also generates a constant which multiplies the resulting function. Thus the 
introduction of a is a logical step. The determination of a will be discussed after the 
SMC equation is derived. 

The closed channel contribution is described by the projected Schrddinger equation, 

(29) 



(1 - aP)H^+) = 0, 

with 

1 

II 


1! 

J* 

1 

and 

1 

KJ 

II 

Note that 

PH = PH a - PV 


= \(HoP + PH 0 ) - PV 
= \(HP + PH)+ l -(VP- 
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Equation (29) can be rewritten as 


[H - HP + PH) + ^{PV - FPpW = 0. (30) 

The SMC equation is obtained by dividing Eq. (30) by a and adding it to Eq. (27), 

j4 (+)^(+) = VSn} ( 31 ) 

with the SMC operator, 

yl(+) = | i(PV + VP) - VG^V +\h- l -(PH + HP)}. (32) 

The SMC equation 4 / uses the LS equation to describe the open channel component 
of 'I'W and the Schrodinger equation for the closed channel component. Because it 
includes both the open and closed channel components, it provides a complete solution 
to the scattering problem. 30 

Using the SMC equation instead of the LS equation, the Schwinger variational 
expression for the T-matrix is 

T mn = (S m |F|*W) + (*L->MS n > - (33) 

If the scattering wave function is expressed in terms of TV + 1-electron trial functions 
/, as in Eqs. (12) and (13), the variational stable expression for T mn is the same as 
Eq* (18), 

T mn = ^(5 m |F|/ r ) J D r .(/ J |V|5„). (34) 

T3 

but with 

(zr 1 ) rj = (/ r |yl<+)|/.). (35) 

3.3. The Projection Parameter a 

The SMC equation is incomplete until a is chosen. A number of ways of determining 
a have been considered in the literature, and they are described below: 

(a) Based on the hermiticity of the principal- value SMC operator. Takat- 
suka and McKoy 18 19 and Lima and McKoy 30 argued that the variational stability of 
T mn requires 

A (+)t = (36) 

In other words, the principal- valued SMC operator A must be Hermitian. Otherwise, 
T mn will be unstable with respect to first order variation of either or 6&-\ While 

Eq. (36) is readily satisfied for trial functions consisting only of L 2 functions, Takatsuka 
and McKoy noted that, if the trial functions included (shielded) spherical Bessel and 
Neumann (or Hankel) functions, Eq. (36) is not valid because the matrix element of 
the kinetic energy operator between Bessel and Neumann functions is non-Hermitian. 
In this case, they showed that if the projection parameter is chosen to be 

a-TV-fl, (37) 

4 / Takatsuka and McKoy’s original derivation 18,19 used the principal- value SMC operator, A , obtained 
by replacing with the principal- value Green’s function, Op. Similarly, the wave functions are 

replaced by those with standing wave boundary conditions. However, most subsequent numerical 
calculations used + h 
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the following matrix element vanishes in the open channel space, 

<P* m \-H-\(PH + HP)\P9 n )=0. 
a l 

The operator \H - \(PH + HP) is the only part of the SMC operator which involves 
the kinetic energy operator. If its matrix element is identically zero in the open channel 
space, Eq. (36) is guaranteed to hold even for a basis which includes spherical Bessel 
and Neumann functions. 

(b) Based on the stability of the T-matrix with respect to first order 
variation of a. Huo and Weatherford 31 pointed out that the SMC equation auto- 
matically incorporates the proper boundary condition through the Green’s function. 
Thus it is unnecessary to include functions with the proper boundary condition in an 
SMC basis. Indeed, almost all SMC calculations carried out so far use an L 2 basis. In 
that case, Eq. (36) is satisfied independent of a. Even if continuum functions are in- 
cluded in the basis, the use of (5-function normalizable continuum functions would again 
ensure the validity of Eq. (36). (A Neumann function is not 6-function normalizable). 
Instead, they considered the role of a as a weight factor for the relative contribution 
between the open and closed channel space. Thus a should be treated as a variation 
parameter and its optimal choice should be based on the variational stability, 

dT mn /da = 0. 

If 'I'Gd and are expanded in terms of Eqs. (12) and (13), the above condition results 
in the following expression, 

EEM/r \H\f.)c.=0. (38) 

r s 

Equation (38) is to be solved together with Eqs. (15) and (16) so that a can be deter- 
mined simultaneously with the expansion coefficients b T and c,. 

Using Eqs. (12) and (13), Eq. (38) can be rewritten as 

(^ n +) \H\^) = 0 . 

Thus it will be automatically satisfied if and/or ^ 1 satisfy the Schrodinger equa- 
tion, 

H 1^ +) = o, 

or 

= 0 . 
m 

Under these circumstances, T mn will be stable for any finite value of a, including 
a = N + 1. This result is related to how the LS and Schrodinger equations are com- 
bined in the SMC equation. A true solution should satisfy both equations, with the 
result independent of how the SMC operator is partitioned. However, in practical cal- 
culations we search for a variationally stable T mn . Equation (38), which is identical to 
the configuration-interaction (Cl) equation in electronic structural calculations, is much 
more amenable to practical calculations than the Schrodinger equation itself. Generally 
we have found that an iterative solution of a, coupled with the calculation of T mn , adds 
approximately 10% to the cost of an SMC calculation with fixed a. 
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(c) Based on supplementing the projected LS equation. An alternate 
proof for Eq. (37) was provided by Winstead and McKoy. 32 Instead of using the SMC 
equation, they looked for the variational stability of T mn by using the projected LS 
equation, Eq. (26), alone. 

T mn = (S m |V|< +) ) + - <*t ) | VP ~ VG^VIfW). 

While the above expression is not variationally stable, they found that stability can be 
achieved by adding to the projected LS equation a term QH. 

T mn = (S m |V|*< + >) + S n ) - - VG<- P +) V + QH\ (39) 

with 

Q = P - R. 

The operator R , appbed to an antisymmetrized N 4 1-electron function, removes the 
antisymmetrization between an N-electron function and the one-electron function for 
the ( N 4- 1 )t h electron. Thus we have 

It can be readily shown that Eq. (39) leads to a variational stable expression for T mn , 
with Eq. (37) for a. However, unlike Eq. (28), the operators P and Q do not span the 
complete space. In view of the fact that P is an iWelectron operator whereas R, and 
hence Q , is an N 4 1-electron operator, we have 

P + Q/l. 

In contrast with the result of Lima and McKoy 30 who used Eq. (28) to prove the com- 
pleteness of the solutions of the SMC equation, the present approach fails to demon- 
strate that the solution of the operator equation 

(VP - VG^V + QH)V m = S, 

is complete. 

3.4. Implementation of the SMC Method 

The trial wave function in the SMC method is expanded by 

*! +) = ££ &£U{<mi • • • N) Xj (N + 1)} + £ 4 n) ©/(i ■■■* + 1). ( 4 °) 

i=i i 1 

where the index i sums over open channel target functions, j sums over the one-electron 
basis set \ used to expand the continuum electron function, and l sums over 0, the 
antisymmetrized N 4 1-electron configuration-state- function (CSF). The expansion co- 
efficients b\^ and b\ n) are to be determined variationally. 

So far, all SMC calculations on molecules have used target functions represented by 
a Cartesian Gaussian basis. Also, they make use of the fact that the correct boundary 
conditions have been automatically incorporated in the SMC equation and employ an 
L 2 basis of Cartesian Gaussian functions to represent the continuum electron. Thus the 
trial form of contains no information of its behavior at the boundary. Instead, 

its initial condition is given by S„ and the outgoing wave by in the SMC 
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equation. For collisions with long range potentials, where an L 2 basis is inadequate, the 
higher partial wave contributions are determined using a Born closure approximation. 41 
In addition, SMC calculations use incoming plane waves, instead of angular momentum 
waves, for the homogeneous solution S . 

S n = $ n (l . . . N)exp(ik n • r N+1 ) 

Thus this is the only method described in this book which uses a linear momentum 
instead of an angular momentum representation for the T-matrix. A plane wave repre- 
sentation is chosen because the two-electron integral between three Gaussian functions 
and one plane wave can be evaluated analytically. 13 

Below, we describe the overall organization of an SMC code, the steps for calculat- 
ing a body-frame fixed-nuclei T-matrix, T mny the transformation from the body-frame 
to laboratory-frame, and the calculation of differential cross sections (dcs). Some of the 
important features in implementing the SMC method are described in more detail in 
subsequent sections. 

(a) Angular quadrature for k m and k n . In the linear momentum represen- 
tation, T mn is a function of both the magnitude and direction of k m and k n , T mn = 
T(k rr ,,k n ). The total energy determines the magnitude of the wave vector, see Eq. (9). 
To simulate the random orientation of a molecule in a gas phase collision, SMC calcula- 
tions for gas phase c- molecule collisions are performed over angular quadratures of k m 
and k n . Thus, instead of positing an electron beam with a fixed direction colliding with 
randomly oriented molecules, we fix a molecule in space and describe its collision with 
electrons coming from different directions, (#,^>). Note that molecular symmetry can be 
employed to reduce the number of quadrature points used. For example, the cylindrical 
symmetry of a diatomic molecule allows us to to calculate the plane wave contribution 
from one tp and deduce the contributions from other <^’s by rotation. Also, part (f) 
shows that the number of partial waves retrieved in a partial wave decomposition of 
T(k m ,k n ) depends on the size of quadrature used. 

(b) Gaussian basis set. Gaussian basis sets for the calculation of molecular 
wave functions have been well studied. A variety of basis sets are available, with vary- 
ing degree of accuracy. All SMC calculations carried out have used bases constructed 
with the segmented contraction scheme, including Dunning’s earlier contracted basis 34 
and his correlated consistent contracted basis. 35 For diatomic molecules and polyatomic 
hydrides, uncontracted bases have also been used. The Gaussian basis for the contin- 
uum electron is obtained by augmenting the target basis with even-tempered diffuse 
functions. Their exponents are determined by 


where ( t r is the exponent of the i th Gaussian of symmetry type T and a r is a constant, 
usually between 2 to 3. The initial (f is taken from the most diffuse function of this 
symmetry type in the target basis. Note that as more diffuse functions are introduced, 
the basis set is closer to redundancy. Hence, instead of placing diffuse functions at 
each nuclear center, it becomes advisable to place the most diffuse functions only at 
the center of mass. Also, if a basis set is truly redundant, the lowest eigenvalue of its 
overlap matrix is zero. This feature can be used as a quick test for basis set redundancy. 

\\ hen the Green s function matrix elements are calculated using the insertion tech- 
nique, an additional basis is needed for the insertion calculation, The choice of this basis, 
which is critical to the success of the insertion method, will be discussed in Sec. 3.5. 
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(c) Open and Closed channel functions. The open channel projection op- 
erator P is defined to span all energetically accessible target functions. However, in 
practical calculations it is frequently not feasible to include all open channels defined in 
this manner, since the size of the calculation will become too large. This problem occurs 
when the electron energy is close to or larger than the first ionization potential of the 
target and the full set of Rydberg states, the ionization continuum, and dissociative 
states become open. Thus a selection process must be made. Generally, a multichannel 
study is undertaken because a particular set of elastic/excitation processes is of interest. 
The open channel space should include the initial and final channels for these processes, 
as well as other open channels which have significant couplings with the channels under 
study. For example, in a study of valence excitations, Rydberg states are generally ne- 
glected because the Rydberg- valence coupling tends to be small. 36 A significant amount 
of experimentation is required in the selection. 

The full Gaussian basis is used to span the continuum electron orbital. This guar- 
antees that the open channel function is invariant when we transform the orbitals to 
achieve an optimal representation of the closed channel space. Notice that, in the SMC 
method, terms of the type given in Eqs. (23) and (24), so called penetration terms, 
belong to the open channel space because they are associated with the open channel 
target functions and the projection operator P does not annihilate them. In the Kohn 
method, which employs the Feshbach projector formalism to partition the open and 
closed channel space, these terms are grouped into the closed channel space instead. 

The closed channel functions are chosen so they can describe, in as compact an 
expansion as possible, the effects of energetically inaccessible channels and/or transient 
negative ions on the scattering process. For a nonresonant process, the closed chan- 
nel space contributes to the description of the mutual polarization effects between the 
target and electron. In the type of trial function used in Eq. (40), this contribution is 
represented by the use of l-hole-2-particle configurations obtained by multiplying the 
dip ole- allowed, singly-excited configurations generated from the target function with 
the continuum orbital. The description “dipole allowed” means the transition dipole 
moment between the singly excited iV-electron CSF and the target function does not 
vanish. The polarization effect will be accounted for if a complete set of such configu- 
rations is included in the trial function. 

In practice, however, this will make the calculation quite large, especially if a large 
one-electron basis is used. Thus truncation of the CSF expansion is necessary. This is 
achieved by employing natural orbitals derived from bound state calculations. Earlier 
SMC calculations 37 used natural orbitals from bound state negative ion Cl calculations. 
Because these are bound state calculations, the continuum electron is simulated by a 
high-lying Rydberg orbital. Nevertheless, the natural orbitals from the Cl calculations 
can be used to provide a shorter expansion for the closed channel functions. A more 
efficient expansion has recently been proposed by Lengsfield et al. using polarized 
orbitals, an approach adopted in recent SMC calculations. 

For processes involving shape resonances, natural orbitals deduced from a bound 
state N + 1-electron Cl still provides an optimal representation of the closed channel 
space. In the cases we have studied so far, the Cl wave function suitable for represent- 
ing the transient negative ion can be readily identified by the presence, in the set of 
natural orbitals for this state, of an orbital with an occupation number of nearly one 
and with a strong antibonding valence character. Based on the size of the correspond- 
ing Cl coefficients, a truncated CSF list can be generated and then used to generate 
the closed channel configurations. This approach has been employed to obtain an ac- 
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curate description of the 2 U g resonance in the elastic scattering of N 2 , 39 and the shape 
resonances in the Be atom and small Be clusters. 22 

Polarization effects are also important for resonant channels so one should generally 
use a combination of the above two strategies. One can either generate a set of polarized 
orbitals from the natural orbitals excluded from the description of the transient negative 
ion, or order the set of discarded CSF’s by the transition dipole moment between the 
AT-electron function associated with the CSF and the target state. A recent calculation 
on elastic scattering of CF 4 , incorporating polarization effects, is based on the second 
approach. Some of the results are presented at the end of part (f). 

The SMC method can also be used to study Feshbach resonances which arise from 
the interaction of the open channel with a bound negative ion associated with an excited 
target state. The closed channel configurations used to represent a Feshbach resonance 
can be generated readily using the bound state Cl techniques described above. An 
SMC study of a Feshbach resonance in e-H 2 collisions has been reported. 40 However, 
this calculation used a frozen core approximation for the excited state orbitals instead 
of Cl natural orbitals. 

(d) Integrals. Three type of integrals are involved: (1) The matrix elements 

(9&\B\9M) and These matrix elements involve integrals between 

Gaussians which are evaluated using existing quantum chemistry packages. (2) The 
matrix elements <¥<“ ) |V P |5 n ) and (S m |F|$W). If V is the nuclear attraction potential, 
it involves integrals between a Gaussian and a plane wave. If V is the two-electron 
repulsion potential, the integrals are between three Gaussians and a plane wave. Ana- 
lytical expressions for both one and two-electron integrals have been derived. 13 42 These 
integrals are complex and their evaluation are more costly than the Gaussian integrals. 
For example, the two-electron integrals are an order of magnitude slower than the cor- 
responding Gaussian integrals. Nevertheless, integral packages are available for their 
calculation. Note that these integrals are calculated over the angular quadrature of 

k n . (3) 1 he Green’s function matrix element. This type of integrals is considered to 
be the bottleneck in a Schwinger calculation. They will be discussed separately in Sec 
3.5. 

(e) Formation of the matrix and its inversion. Once the integrals are 

calculated and the open and closed channel configurations are chosen, the gathering 
of the matrix elements follows the structure of quantum chemistry codes. The 

operation of the projection operator P is also straightforward. For SCF target functions, 
this operation is done directly but for Cl target functions, it can be more efficiently 
performed using density matrices. In Sec. 3.6 we shall consider how SMC calculations 
avoid a certain type of pseudoresonance when correlated target functions are used. 

Due to the use of complex boundary conditions, the A ^ matrix elements are 
complex. Its inversion is done using matrix inversion routines for complex matrices 
available in many computational science libraries. The body-frame T(k m ,k n ) is then 
obtained from Eq. (34) using the inverted A matrix, If the projection parameter a 
is to be determined variation ally, then an iterative solution of a is coupled with the 
calculation of T(k m ,k n ) at this step. 

(f) Frame transformation and cross section expression. To obtain the 
differential cross section, we need to transform the body-frame T-matrix into the 
laboratory-frame. Let ( 8 n , p n ) and (0 m , <p m ) be the angular coordinates of k n and 

in the body-fixed frame. We first expand T(k m , k n ) in terms of partial waves of (0 m , 
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(41) 


V^m ) 

T(k m ,k n ) = k n , 

The body-frame partial wave T-matrix is determined by 

0 n ) = J d9 rn siu9 rn J d<p m Y lfi (9 m , y? m )T(k m , k n ). 

Next the body-frame partial wave T-matrix is transformed to a laboratory-frame with 
the new Z-axis in the direction of k n . The Euler angles for the frame transformation is 

(0 ,0 n ,v? n ). 

T L (9, p, 9 n , (/?„) = Tln(k m , A: n , (/? n > V^n)- 

Here (0,<p) is the solid scattering angle in the laboratory frame. Then |T i | 2 is averaged 
over 9 n ,p n to account for the random orientation of the target, 

a(6,(j>) — f d6 n sin6 n f d<p n \T L {6,ip,<p n ,8n)\ 2 ■ 

lOTT * 5 K n Jo Jo 

The physical cross section is obtained by averaging over the azimuthal angle p, 

Integration of cr(9) over 9 gives the integral cross section. The integral cross section can 
also be obtained directly in the body-frame, 

hv I dkn I 

It should be mentioned that the number of partial waves deducible from Eq. (41) 
is determined by the size of the angular quadrature used to calculate the body-frame 
T-matrix for various plane wave orientations. While in principle each plane wave can 
be decomposed into an infinite sum of angular momentum waves, the finite angular 
quadrature used in the calculation limits the number of angular momentum waves 
obtainable from such decomposition. However, as discussed in part (g), the scattering 
of high l partial waves can be described by perturbation theory and does not require 
a full scale variational calculation. Thus the limited number of partial waves deduced 
from Eq. (41) does not present a problem. In most cases, a (8x8) or (10x10) Gaussian- 
Legendre quadrature for k m and k n give l max = 6 for the angular momentum wave. 

A recent calculation of e-CF 4 elastic scattering, 43 including polarization effects, 
serves as an example of an SMC calculation which employs a set of optimized closed 
channel configurations. An elastic scattering calculation in the static-exchange (SE) 
approximation has been reported by Winstead et al. 44 who improved an earlier SE cal- 
culation by Huo 45 by employing a larger Gaussian basis. They reported an A\ resonance 
at % 13 eV and a T 2 resonance at « 11 eV. The present polarization calculation used 
a Gaussian basis of 10s6pld functions at each nuclear center. The target function was 
described by an SCF function at the experimental equilibrium geometry and an (8x8) 
angular quadrature for the wave vector was used. The closed channel space was chosen 
using the procedure described in part (c) above. Thus bound state Cl calculations were 
carried out for negative ions of 2 A U 2 T 2x , 2 T 2y and 2 T 2z symmetries and the Cl wave 
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Electron Energy (eV) 


Figure 1 . Eigenphase sums for electron-CF 4 scattering, showing resonances in the 2 A\ (dashed line) 
and 2 T 2 x (solid line) partial channels. 



Figure 2. Differential cross section for elastic electron-CF 4 scattering at 8 eV. Solid line, present po- 
larized result, filled circles, experiment of Mann and Linder, 46 and open circles, experiment of Boesten 
et al. 47 



Scattering Angle (Degrees) 


Figure 3. Differential cross section for elastic electron-CF 4 scattering at 11 eV. Solid line, present 
polarized result, and filled circles, experiment of Mann and Linder. 46 
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functions with a strong antibonding character were identified as good candidates to 
describe the transient negative ions. It was found that the natural orbitals from the 2 A\ 
ion calculation could also provide rather compact representations of the Cl functions 
of the 2 T 2z , 2 T 2y and 2 T 2z ions. Thus the 2 A l natural orbitals were used in the scat- 
tering calculation and lists of CSF's for the four symmetries were generated from the 
respective negative ion Cl calculations. The lists were augmented to include polarization 
effects based on the size of the transition moment. They were then used to represent the 
closed channel functions in the polarized scattering calculation. Figure 1 presents the 
eigenphase sum of the A x and T 2x partial channels. The position of the A\ resonance is 
lowered to 9.0 eV and the T 2 resonance to 8.6 eV, in much better agreement with ex- 
perimentally observed features. 46 47 Figures 2 and 3 present the theoretical differential 
cross section (dcs) at 8 and 11 eV, together with experimental data. 

(g) Born closure for long range potentials. When the interaction potential 
or transition potential includes long range multiple potentials, such as the dipole poten- 
tial in elastic scattering or the dipole transition potential in dipole-allowed transitions, 
the number of partial waves required to converge the cross section, especially in the for- 
ward direction, is significantly larger than what can be provided by an SMC calculation 
with moderate size angular quadrature. On the other hand, due to the large centrifugal 
barrier for the higher partial waves, they are prevented from reaching the inner region 
of the target and the scattering tends to be weaker. Thus the Born approximation is 
expected to be valid for the high l partial waves and can be employed to describe their 
contributions 48 50 to the cross section. Gibson et al. 41 made use of the following Born 
closure approximation, 

<t(0) = a FBA {6) + [<t smc (9) - a F F B E A {e)}. (42) 

to include the high l contributions in the dcs. Here F B A denotes the first Born approx- 
imation and <j£f A (0) is obtained from a finite expansion of the first Born cross section 
containing exactly the same number of partial waves as a SMC (0). In Eq. (42), the Born 
value for the low l contributions to the dcs is subtracted out. Thus we use the SMC 
result for the low l and Born value for the high l contributions to the dcs. In evaluating 
a FE A (Q)i care should be taken so that the same angular quadrature of k m and k n are 
used as in the SMC calculation. 



Figure 4. Differential cross section for A rl E+ -» transition of H 2 at 15 eV. Solid line, SMC’ 

result, dashed line, distorted wave result, 48 and crosses, experimental data of Srivastava and Jensen. 49 
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Figure 5. Differential cross section for — ► B 1 E+ transition of H 2 at 20 eV. Solid line, SMC 

result, dashed line, distorted wave result, 48 crosses, experimental data of Srivastava and Jensen, 49 and 
triangles, experimental data of Khakoo and Trajmar. 50 


Equation (42) has been applied to a two-channel SMC calculation to study the 
dipole allowed transition — > R 1 !!* in H 2 . Figures 4 and 5 present the 15 and 20 

eV dcs from the SMC calculation with Born closure, the distorted wave calculation of 
Rescigno et al., 51 and the experimental data of Srivastava and Jensen 52 and Khakoo 
and 1 raj mar. 

It should be pointed out that in their implementation of the Born closure approxi- 
mation to study elastic scattering of polar molecules, Rescigno et al. 4 54 took a slightly 
different approach and applied it to the body-frame T-matrix instead of the dcs. 

3.5. Evaluation of the Green’s Function Matrix Elements 

The Green’s function matrix element 

M = (¥J\VG { + ] V\¥ n ^), 


can be expressed as 


M = Mr + iMj. 

The imaginary part Mi is given by the residue, 

Mi = ktf 


(43) 


where l sums over the poles of Gp^. Since G p^ is a projected Green’s function, this is 
equivalent to a sum over the open channels and ki is just the wave vector associated with 
the open channel 1. Also, (5/| V r |^ r { l + ^) is the same matrix element used in the expression 
of the body-frame T-matrix, T\ n . The integration over kj in Eq. (43) can use the same 
quadrature as the body-frame T-matrix. Consequently, the calculation of Mi involves 
very little extra work. 

The real part of M is given by the principal- value integral, 


Mr = 


-J—V [ °°de E 

( 2*) 3 Jo ^JE 


klL 


Si) — e -f iS 


(44) 
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with 


u mn (k) = I (45) 

Here e = \k 2 . Due to the integration over e and k, Mr is much more difficult to 
compute. Its efficient and accurate calculation is an important factor in the development 
of a practical SMC method. Two approaches to computing this integral are described 
below. 


(a) Direct numerical calculation over a quadrature of k and e. An opti 
mized quadrature is obviously important in this approach. To avoid the rapid variation 
of the integrand in the e coordinate near the poles of the Green’s function, we follow 
the work of Walters 55 and Heller and Reinhardt 56 and separate the e-integration into 
two regions, with the poles of the Green’s function all located in the first region. The 
rapid variation in the first integrand is avoided by the addition and subtraction of a 
term such that the numerical integration can be carried out over a smooth integrand. 
The e-integral in the first region is rewritten as, 

_ f €a J k Um n (k) f ta j f kU mn (k) klUmn{kl) k{ Umn{kl) | 

(46) 

with e a > (E — Ei) for all l. In Eq. (46), the second and third term in the integrand are 
obtained from the first by replacing k and k with ki and ki . Because U m n(ki) and ki are 
independent of e, the last integral on the right-hand side of Eq. (46) can be evaluated 
analytically and Eq. (46) is rewritten as, 


h 


£ {kiU^k^lni 

l£opcn 


^ a 



de 


kU mn (k) 

E - E t - e 


k Umnjkl) 1 1 

E-Ei-e J J 


Since the first and second term in the e-integand above have poles at the same place, the 
large changes in the integrand near the pole region are canceled out and the integrand 
remains smooth, suitable for numerical integration. 

To demonstrate the convergence of the numerical quadrature scheme, we calcu- 
lated the elastic scattering of CO in the static-exchange approximation, using an SCF 
target at the experimental equilibrium geometry. A Gaussian basis, with 5s3pld con- 
tracted functions at each nuclear center, represented the continuum electron. The SMC 
calculation included up to l max = 6. Born closure was not included since our purpose 
was to demonstrate the convergence of the VG^V quadrature, and the introduction 
of a constant correction term from the Born closure would not affect our conclusions. 
Also, all calculations used the Gauss-Legendre quadrature. 

Figure 6 presents the partial integral cross section for e-CO elastic scattering. A 
(20x20) quadrature was used for the angular integration and a 20-point quadrature for 
the second region e integral. Three different quadratures were used for the first e region, 
20, 32, and 48 points. The three cross section curves virtually coalesce, indicating a well 
converged result. 

The cross section is much more sensitive to the quadrature used in the second € 
integration. Figure 7 shows the cross sections obtained using 20, 32, and 48 quadrature 
points for the second region integration and keeping the other quadratures constant. 
Note that the peak of the cross section curve at fs 3.5 eV, arising from the 2 II shape 
resonance in CO, is shifted to a lower energy when a larger quadrature is used. Also, the 
32- and 48-point calculations have similar resonance positions, but the cross sections 
themselves differ by ~ 10% at the peak. 
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Sensitivity to the angular quadrature was also studied. Figure 8 presents the results 
for the variation in the 0 quadrature, keeping all other integration parameters constant. 
While the 32- and 48-point integrations give substantially similar results, the 20-point 
quadrature has not yet converged. In particular, note that the maximum of the cross 
section, at % 3.5 eV energy, is shifted to a higher energy when a larger quadrature is 
used. 



Energy (eV) 


Figure 6. Partial integral cross section (/ mai =6) for e-CO elastic scattering in the static-exchange 
approximation. A (20x20) quadrature is used for the angular integration in the VG^V term and 
20-point for the second region e-integration. The quadrature size for the first region e-integration is: 
circle, 20 points, square, 32 points, and triangle, 48 points. 



00 1.0 Z.O 3 0 4.0 SO 


Energy (eV) 

Figure 7. Same as Fig. 6. A (20x20) quadrature is used for the angular integration and 20-point 
for the first region e-integration. The quadrature size for the second region e-integration is: circle, 20 
points, square, 32 points, and triangle, 48 points. 

When a bigger basis set is used for the target function, the demand for a larger 
angular quadrature becomes higher, particularly if the basis set is uncontracted and 
has very large exponents. Also, the second e integral shows more sensitivity to the size 
of the angular quadrature. These features suggest that the off-shell contribution to the 
VG^pW term, coming from the high energy region of the target spectrum, may be 
responsible for the need of large quadratures. In the above example, a converged result 
requires over forty thousand quadrature points. 
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Figure 8. Same as Fig. 6. A 20-point quadrature is used for the first region e-integration, 48 for the 
second, and 20 for the (p integration. Circle, 20-point 9 quadrature, square, 32 points, and triangle, 48 
points. 


The direct numerical integration approach was first investigated by Lima et al. 5 ' 
To expedite the integral calculation, Winstead et al. 58 have successfully adapted the 
integral calculation and the three-index transformation 5 ^ to the highly parallel com- 
puters. Parallel computers are particularly well suited for this purpose because each 
integral calculation is independent of others. SMC calculations using this technique 
have been carried out on many polyatomic molecules. 12 

(b ) Insertion using a Gaussian basis. Provided that a Gaussian basis satisfies 
the approximate relation 

E IP«)(Pa| ~ 

a 

the free-particle Green’s function can be projected into this basis, called the insertion 
basis. The principal- value Green’s function matrix element is now written as 

= E EE^imP-tr ^ ))(*, p b (v N+i )\V\¥Y) 

' lEopen a b 


x V 


j™ de J dk k 


( Pa (r^ +1 )le’ k r ^'){e- k r ^M P6( r^ +1 )) _ 
( L — t/) — t f 2^ 


(47) 


Both e and k integrals are expressed in closed form. 59 Furthermore, calculations of these 
Gaussian integrals are an order of magnitude faster than the three-Gaussian-one- plane 
wave integrals in Eq. (45). 

It is important to recognize that in Eq. (47), an L 2 basis is used in the evaluation 
of an integral over the Green’s function, but not to represent the Green’s function itself. 
Just as an L 2 basis can be used to represent a trial function in a Schwinger variational 
calculation because the trial function is always multiplied by V in the expression of 
the 7-matrix, in the present case the free-particle Green’s function is integrated over 
w ith the trial \fd+) itself expressed in terms of another L 2 basis. Thus, the 
long range oscillatory part of the Green’s function does not contribute to the integral. 
In this respect, the application of an L 2 basis to evaluate Mr is analogous in spirit to 


5 ' The fourth index is the plane wave. 
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Schneider’s calculation 60 of the Schwinger separable potential by projection onto an L 2 
basis. 

The insertion technique was used in early SMC calculations when large, uncon- 
tracted Gaussian bases were used for the scattering basis, with the insertion basis 
usually composed of the scattering basis supplemented by additional Gaussians. Good 
convergence behavior was observed. 39 The change to the direct numerical integration 
technique by Lima et al. 57 was motivated by the convergence problems they encoun- 
tered. Such problems appeared to be related to choosing a suitable insertion basis. 

In view of recent progress made in the field of quantum chemistry in developing 
flexible Gaussian bases, 33 and the efficiency in the computation of Gaussian integrals 
themselves, we re-investigated the insertion technique. 20 The choice of a suitable basis 
was guided by the convergence characteristics observed in the numerical integration 
approach, namely that the off-shell high energy contribution to Mr seems to play an 
important role in the convergence. To describe this type of contribution, we included 
Gaussians with large exponents. Also, a multicentered Gaussian basis readily simulates 
the dependence on the angular quadrature. When the insertion basis was chosen using 
these guidelines, good convergence was observed. As an example, we studied the e-CO 
static-exchange problem using the same scattering basis as discussed above. The inser- 
tion basis is obtained by supplementing the scattering basis by 10sl8p9d Gaussians at 
each nuclear center and 16s5p9d Gaussians at the center of mass. 6 / Figure 9 compares 
the cross sections obtained using the insertion method versus the direct numerical in- 
tegration approach. Good agreement is observed. The largest discrepancy between the 
two results, less than 3%, is found at the low energy region. In the 2 II shape resonance 
region, the two curves almost coalesce. 



Energy (eV) 


Figure 9. Partial integral cross section (/ maz — 6) for e-CO elastic scattering in the static-exchange 
approximation. Square, calculated using a direct numerical integration, and circle, calculated with the 
insertion technique. The numerical integration uses a (32x32) angular quadrature, 20-point for the first 
€ integration, and 48-point for the second. The insertion basis is described in the text. 


While the insertion basis used for the CO calculation appears to be very large, it 
should be pointed out that the insertion basis has 238 Gaussian functions and the total 
number of Gaussian integrals used is 15.2 million. The numerical integration calculation 

6y/ Not all components of the d functions were used. Also, a smaller size p* basis was used at the 
nuclear center. 
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requires more than 2.78 billion three-Gaussian-one-plane wave type integrals. Taking 
into account the difference in the cycle time required to calculate the two types of 
integrals, the insertion technique requires orders of magnitude less CPU and storage 
resources. 

Recently, Bettega et al. 61 applied the pseudopotential method and removed the 
core electrons from the SMC calculations. This raises an interesting question. Since the 
off-shell high energy contribution to Mr that we observed originate, at least partially, 
from the core electrons, the use of pseudopotentials may lessen the demand for large 
quadratures in the numerical integration technique and the demand for tight Gaussians 
in the insertion technique. This aspect remains to be tested. 

3.6. Correlated Target Function in the SMC Method 

Because molecular bonding changes with nuclear geometry, multiconfiguration 
wave functions are required to provide a consistent description of the target over a 
range of geometries. Similarly, the description of several electronic states with a com- 
mon set of orbitals, a practice frequently used in the study of electronic excitations 
by electron or photon impact, requires the use of multiconfiguration target functions. 
Thus correlated target functions are used with some frequency in ab initio studies of 
e-molecule collisions. 62 65 

An ubiquitous feature in the use of multiconfiguration target functions is the pres- 
ence of pseudoresonances at intermediate energies. 62, 65 Lengsfield and Rescigno 62 at- 
tributed these resonances to the fact that multiconfiguration target functions can in- 
troduce certain terms in the N + 1-electron wave function which are associated with 
excited states of the target excluded from the open channel space. Pseudoresonances 
appear at an energy where the excited target state becomes open. Consider the example 
of a two-configuration function for the ground state of H 2 8 / 

$x = a\4>\ + 2 , 

with 

4 >\ = -4{l^ 9 a(l) lcr 9/?( 2 )}> 

<t> 2 = .4{l<T u a(l)l<7-„/?(2)}. 

Here is the simplest function that can describe the dissociation of H 2 properly. 

Two N + 1-electron CSF’s are used to relax the orthogonality between the target 
and the continuum electron, 

0a = .4w + ,{*x(l,2)l<r t a(3)} = a*4 w+l {fc(l,2)l«- B a(3)}- 

and 

= *4jy+ij$x(l> 2)l<r u a(3)j = ai*4jv+i 2)l<r u a(3)|. 

Here An+\ antisymmetrizes electron 3 with electrons 1 and 2. Note that © a includes the 
terms \<T g a(l)l<T u a(2)la u (3(3) and la g <x(l)l<j u /3(2)l<r u (x(3). These terms are associated 
with the 6 3 E+ target state, with the la u orbital representing the continuum electron. In 
the Kohn calculation it was found that if the 6 3 £* state was not included in the open 

7 / As a result of symmetry operations, in both cases the actual number of distinct integrals to be 
calculated are less. However, the ratio of the the two numbers will be approximately the same when 
symmetry is taken into account. 

8// This example was originally used by Lengsfield and Rescigno. 
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channel space, a pseudoresonance would occur around the energy where this channel 
became open. Similar pseudoresonances were also observed in the R-matrix calculation 
of the electronic excitations of H 2 . 65 The strategy used in the Kohn method to eliminate 
such terms has been discussed. 11 62 

Due to the use of the iV-electron projection operator, the SMC method treats corre- 
lated target functions differently. As seen in Sec. 3.1, the iV-electron projection operator 
selects the terms that are associated with the target function from an antisymmetrized 
N + 1-electron CSF, resulting in an un-antisymmetrized product of the JV-electron 
target function and an one-electron orbital describing the continuum electron. For a 
multiconfiguration target function, this operation has additional consequences. As an 
illustration, consider the same two-configuration wave function. The projection operator 
P is given by 

|*x>(*x|. 


The operation of P on 0 a gives, 

P|0 a ) = ^|^(l,2)l<7 9 a(3). 

Thus the operation of P recovers the (f> i component in 0 a , originally eliminated by 
antisymmetrization. In this manner, the projector P removes the character from 

0 a , and P0Q, unlike 0 a itself, consists only of the open channel target. 

In an SMC calculation, a pseudoresonance usually appears when the denominator 
in the expression of the T-matrix becomes particularly small at a certain energy. Because 
0 a includes a term that behaves like a open channel function, the matrix element 

as a function of energy E goes through a minimum around the energy when that channel 
becomes open. On the other hand, PG a does not have any 5 3 E + character, and 

does not exhibit the same minimum. 

Based on this analysis, an SMC calculation is expected to perform differently from 
other methods when a multiconfiguration target is used. Figure 10 presents an SMC 
study of the 2 £+ channel elastic scattering of H 2 , using the static-exchange approxima- 
tion. The calculation was carried out at R = 4.0 a 0 , using a Gaussian basis of 6s6p at 
H and 4s4p at the center of mass, The Cl coefficients of the two-configuration target 
function are 0.8785 and -0.4778. As seen in Fig. 10, the SMC cross section curve is 
smooth over the energy range of 0.25 to 21 eV. The pseudoresonance found in other 
methods, arising due to the term 0 Q , is absent here. Thus it appears that the role of 
the projection of P on 0 a is to associate 0 a with the correct open channel, and thereby 
removes the spurious resonance. 

The above results suggests that the jV-electron projection operator has additional 
advantages besides providing a convenient treatment for exchange. The use of correlated 
target function in the SMC method is relatively recent. Only a few studies have been 
made using a multiconfiguration target function. 22 Future studies should explore this 
aspect further. 
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Figure 10. Integral cross section for e-H 2 2 E+ channel elastic scattering in the static-exchange ap- 
proximation, calculated using a twoteim Cl function for the target 


4. THE USE OF SMC RESULTS IN THE STUDY OF VIBRATIONAL 
EXCITATIONS 


All SMC calculations use the fixed-nuclei approximation. However, the fixed-nuclei 
SMC results can be readily used as a starting point in the application of the nuclear 
dynamics treatment of Domcke. Domcke’s treatment, which is based on the projection 
operator formalism, has been reviewed recently, 23 and it will not be repeated. Here 
we consider an example of combining the two approaches in the study of the resonant 
enhanced vibrational excitation of N 2 by electron impact. Because the 2 U g resonance 
in N 2 is a pure d-wave resonance, we rewrite the T-matrix for the Z= 2 partial wave in 
terms of the A"-matrix, 

K(E, R) 


where R is the internuclear distance. The A'-matrix element is related to the width 
T(EyR) and shift A (E,R) functions by 


K(EyR) = - 


£ (E,R) 

2{E - e d (R) - A( A, R)} ' 


Here e d is the difference between the target and negative ion potential curves. The shift 
function is given by the principal value integral, 


*<*•*> - 5 h v l dE 


, r(E',A) 

E - E* ’ 


and the width function can be expressed in terms of the potential UdE which Domcke 
called the entry amplitude, 

T{E,R) = 2it\lJ dE \ 2 . 

An analytical fit of U d E can be obtained using the A-matrix determined from the SMC 
calculation. Notice that U dE is a nonlocal potential, depending on both R and E. 

The calculation on N 2 has been reported previously. 39 66 In applying Domcke’s 
method, we used a Schwinger- type separable potential approach to evaluate the T- 
matrix element. A complete set of the bound vibrational wave function of the transient 
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negative ion, obtained using a numerical solution of the 1-d schrodinger equation, was 
used in the separable potential calculation. The negative ion curve was determined 
empirically by Berman et al. 67 

The results of this calculation compared well with experiment, and the comparison 
will not be repeated. Here we consider an example of applying our result to an experi- 
mentally inaccessible region. In very high temperature plasmas, as in the flow field of a 
high-speed space vehicle upon entry to the earth’s atmosphere, molecules are very hot 
ro- vibration ally and it is difficult to duplicate such conditions in the laboratory. Figure 
11 presents the v = 0 — ► 1 and v = 3 — > 4 vibrational excitation cross sections of N 2 
by electron impact. While many experimental measurements on transitions from v = 0 
have been done, there has been only one reported measurement with the initial state at 
the v = 1 level. 68 Due to the difficulties involved in preparing a significant quantity of 
N 2 at the v — 3 level, no measurement has been reported for excitations from this level 
and theoretical data is the only available source for these cross sections. It is seen from 
Fig. 11 that the resonance structures for the v = 3 — > 4 transition is more extended 
both at high and low energies in comparison with the v= 0 — ► 1 transition, and not as 
strongly peaked. The calculation assumed the molecule was initially at J = 50, with the 
centrifugal barrier in the vibrational potential determined by this J value. Note that 
our calculation did not include rotation in the dynamics and the effect of rotation only 
entered through the centrifugal term. 



Electron Enerty (eV) 


Figure 11. N 2 vibrational excitation cross sections by electron impact. Solid line, v = 0 — » 1 , dashed, 
v = 3 — > 4. The calculation assumes the molecule is at J = 50. 


Figure 12 presents the v = 0 — + 1 excitation cross sections for N 2 at J= 0, 50, and 
150. Again rotational effects contribute only through the centrifugal barrier of vibra- 
tional motion. Because the differences between the potential curves of the neutral and 
transient negative ion, the large centrifugal barriers at J — 50 and 150 have differential 
effects on the vibrational excitation. The resonance structure in the cross section moves 
to lower energies and narrows as J increases. 

Future studies that couple vibrational and rotational motion at high J will extend 
our understanding of the collisional processes involving ro-vibrationally hot molecules. 
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Figure 12. N 2 v = 0 1 vibrational excitation cross section by electron impact. Solid line, J — 0, 

dotted, J — 50, and dashed, J= 150. 

5. SUMMARY 

The SMC method, a modified version of the SV method, has proven to be a ver- 
satile tool in the study of e-molecule collisions. The SMC method has been successfully 
applied to the study of technologically important electron collision processes involving 
polyatomic targets. This application will continue as the SMC method is further devel- 
oped and refined. The present discussion of its implementation also serves to illustrate 
some of the physical basis in its formulation, and provides a foundation for the future 
development of the method. 
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